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We perform a renormalisation group (RG) analysis of neutrino masses and mixing 
angles in the see-saw mechanism in the minimal supersymmetric standard model with 



three right-handed neutrinos, including the effects of the heavy neutrino thresholds 



We focus on the case that one of the right-handed neutrinos provides the dominant 
contribution to the 23 block of the light Majorana matrix, causing its determinant to 
approximately vanish and giving an automatic neutrino mass hierarchy, so called sin- 
gle right-handed neutrino dominance which may arise from a U(l) family symmetry. 
In these models radiative corrections can increase atomospheric and solar neutrino 
mixing by up to about 10% and 5%, respectively, and may help to achieve bi-maximal 
mixing. Significantly we find that the radiative corrections over the heavy neutrino 
threshold region are at least as important as those usually considered from the lightest 
right-handed neutrino down to low energies. 
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1 Introduction 



The latest atmospheric neutrino results based on 1117 days of data from Super 
Kamiokande are still consistent with a standard two neutrino oscillation — > v T 
with a near maximal mixing angle sin 2 2#23 > 0.88 and a mass square splitting Am^ 
from 1.5 x 1CT 3 to 5 x 1CT 3 eV 2 at 90% CL [I]]. The sterile neutrino oscillation 
hypothesis — > v s is excluded at 99% CL. Super Kamiokande is also beginning to 
provide important clues concerning the correct solution to the solar neutrino prob- 
lem. The latest results from 1117 days of data from Super Kamiokande sees a one 
sigma day-night asymmetry, and a flat energy spectrum, which together disfavour the 
small mixing angle (SMA) MSW solution ||, the just-so vacuum oscillation hypoth- 
esis 0] and the sterile neutrino hypotheses. All three possibilities are now excluded 
at 95% CL. The results allow much of the large mixing angle (LMA) MSW region, 
which now looks like the leading candidate for the solution to the solar neutrino prob- 
lem. For example a typical point in the LMA MSW region is sin 2 29 12 ~ 0.75 and 
Am 2 2 « 2.5 x 1(T 5 eV 2 f§. 

The see-saw mechanism || implies that the three light neutrino masses arise from 
some large mass scales corresponding to the Majorana masses of some heavy "right- 
handed neutrinos" N R M R 9 R (p,q = 1, • • • , Z) whose entries take values which extend 
from ~ 10 14 GeV down to perhaps several orders of magnitude lower. The presence of 
electroweak scale Dirac mass terms rnl R (a 3 x Z matrix) connecting the left-handed 
neutrinos v\ (i = 1, . . .3) to the right-handed neutrinos N R then results in a very 
light see-saw suppressed effective 3x3 Majorana mass matrix 

m LL = m LR M RR m T LR (1) 

for the left-handed neutrinos u l L , which are the light physical degrees of freedom 
observed by experiment. If the neutrino masses arise from the see-saw mechanism 
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then it is natural to assume the existence of a physical neutrino mass hierarchy 
m Vl <C m„ 2 <C Tn U3 , which implies Am| 3 ~ m l 3 i an d Amf 2 ~ m ^ 2 ' wm ch fixes 
m U3 « 5.6 x 10~ 2 e\^, and (assuming the LMA MSW solution) m^ 2 w 5.2 x 10~ 3 el^, 
with rather large errors. Thus m U2 /m U3 ~ 0.1. In view of such a 23 mass hierarchy 
the presence of a large 23 mixing angle looks a bit surprising at first sight, especially 
given our experience with small quark mixing angles. Several explanations have been 
proposed 0, but the simplest idea is that the contributions to the 23 block of the light 
effective Majorana matrix come predominantly from a single right-handed neutrino, 
which causes the 23 sub determinant to approximately vanish. This mechanism, called 
single right-handed neutrino dominance (SRHND), was proposed in ||, and developed 
for bi-maximal mixing in ||. 

In this paper we shall be concerned with the effect of radiative corrections in 
models based on SRHND and U(l) family symmetry ||, and in particular on the 
case of bi-maximal mixing in these models |§ . We choose these models because 
analytic estimates suggest that the 23 hierarchy arises from a physical mechanism 
which accounts for the smallness of the 23 subdeterminant, and so this hierarchy 
should be stable under radiative corrections. Although the radiative corrections to 
atmospheric mixing are only a modest 10%, this can nevertheless play an important 
role in achieving near maximal atmospheric mixing. Our results may be compared 
to the several RG studies of various models which already exist in the literature [ |T0"|| , 



TT|. Many of the existing studies do not take proper account of the heavy right- 
handed neutrino mass thresholds, often only including the effects of running due to 
the low energy dimension 5 see-saw operator from the see-saw scale down to low 



energies |T(| , or running from the high energy scale down to low energies including a 



single right-handed neutrino mass threshold |TTJ, assuming that all the right-handed 
neutrinos are degenerate. Our results indicate that the corrections in running from 
the high energy scale, through the heavy neutrino threshold region, down to the 
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lightest right-handed neutrino mass scale are just as important (and in some cases 
more so) than the traditionally calculated radiative corrections in running from the 
lightest right-handed neutrino mass scale down to low energies. 

The paper is organised as follows. In section 2 we define the MSSM with Z right- 
handed neutrinos, where the heavy neutrino mass matrix arises from the vacuum 
expectation value (VEV) of a singlet field S, and write down the renormalisation 
group equations (RGEs) relevant for running from the unification scale Mjj ~ 2 x 
10 ie GeV down to low energies. In section 3 we review the analytic discussion of 
the phenomenological conditions for SRHND and LMA MSW involving the different 
types of heavy Majorana neutrino textures, and show how a £7(1) family symmetry 
may be used to satisfy them [g]. In section 4 we discuss three explicit examples of this 
kind, and then perform our numerical RG analysis of these cases. Section 5 concludes 
the paper. 

2 The MSSM with Z Right-handed neutrinos 

We consider the Yukawa terms with two Higgs doublets augmented by Z right-handed 
neutrinos, which, are given by 

C yuk = e ab [-Y^ H^Q^Uj + Y i d j H° l Q b i D c j + Y^H^L^Ej - Y i u p H"L b i Np 

+ 1 -y^n;n^\ + h.c. (2) 

where e a b = —£ba, ei2 = 1, and the remaining notation is standard except that the 
Z right-handed neutrinos have been replaced by their CP conjugates N£ with 
p, q = 1, . . . , Z and we have introduced a singlet field S whose vacuum expectation 
value (VEV) induces a heavy Majorana matrix Mrr =< S > Y&r. When the two 
Higgs doublets get their VEVS < H\ >= v 2 , < H\ >= v\ with tan/3 = v^/v\ we find 
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the terms 

C yuk = vzYVUiU? + v l Y l d j D i D c j + v 1 Y^ j E i E] + v 2 Y» p N t N c p + ^M p R \N^N c q + He. (3) 

Replacing CP conjugate fields we can write in a matrix notation 

C yuk = U L v 2 Y u U R +D L v 1 Y d D R +E L v 1 Y e E R +N L v 2 Y v N R +^M RR N R +H.c. (4) 

where we have assumed that all the masses and Yukawa couplings are real and written 
Y* = Y. The diagonal mass matrices are given by the following unitary transforma- 
tions 

vzYdiag = VuLV 2 Y u Vl R = diag(m u , m c , m t ), 
ViY£ ag = V dL v x Y d Vl R = diag(m d , m s , m b ), 
v i Y dia 9 = V eL v x Y e V} R = diag(m e ,m M ,m r ), 
M RR 9 = n RR M RR n RR = diag(M R1 , . . . , M RZ ). (5) 

Below the mass of the lightest right-handed neutrino M R ± the right-handed neu- 
trino masses may be integrated out of the theory, which corresponds to replacing the 
last two terms in Eq.^| by a dimension 5 operator 

- e^H^N; + \m™ r N;NI + H.c. -> - l -K l3 {t ab H a u L\){e al yH^L b ;) + H.c. (6) 

where k = Y u M RR Yj is simply related to the see-saw mass matrix in Eq.|l| when the 
Higgs fields are replaced by their VEVs 

m LL = v\k. (7) 

Having constructed the light Majorana mass matrix it must then be diagonalised by 
unitary transformations, 

m LL 9 = VuLm LL Vl L = diag(m^, m^, m V3 ). (8) 
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The CKM matrix is given by 

Vckm = V uL V] L (9) 



The leptonic analogue of the CKM matrix is the MNS matrix defined as [12 



Vmns = V eL Vj L . (10) 

which may be parametrised by a sequence of three rotations about the 1,2 and 3 axes, 
as in the standard CKM parametrisation, 

Vmns — -R23-R13-R12 (11) 



where 



R 



23 



/ 1 \ / C13 S13 

c 23 s 23 , R 13 = 1 I . /? l2 

V -S 23 C 23 / V -Sl3 C13 



/ C12 S12 

-S12 Cia I (12) 

V 001 



where s^- = sin 6^, Cij = cos%, and 6^- refer to lepton mixing angles. Note that we 
completely ignore CP violating phases in this paper. From the unitarity conditions of 
the MNS matrix elements in Eq.ll, and parametrisation in Eq.12, the mixing angles 
can be expressed in terms of elements of Vmns as 

AV 2 V 2 

S soi = sin 2 20i2 = (y2^y\y > (I 3 ) 
4V\V\ 

s "= sin2 ^=wf^ (14) 

In principle these mixing angles may not be exactly the same as those obtained from 
two generation analysis of the experimental results. 

The renormalisation group equations (RGEs) to one-loop order are: 
dY u 1 



dt 16tt 2 
dY d 1 



dt 16vr 2 
dY v 1 



dt 16tp 



[N g .Y u + Y u .N u + (N Hu )Y u ] 
[N q .Y d + Y d .N d + (N Hd )Y d ] 
[N l .Y u + Y".N v +(N Bu )Y v ] 



dY e 1 

= -[N l .Y e + Y e .N e + (N Hd )Y e ] 



dt 167T 

dY RR 1 



■[N u .Y RR + Y RR .N u + {Nx)Y RR ] 



dt 167T : 



where the wavefunction anomalous dimensions are 



N q = {lgl + \gl + ^gl)i-Y u Y^-Y d Y d1 

N u = {\gl + ^gl)I-2Y^Y u 

N d = (lgS + lg*)I-2Y*Y* 

Ni = {\gl + ^gl)I-Y«Y*-Y»Y^ 

N e = £g 2 1 )I-2Y^Y e 
5 

n v = -2r f r - y^y rr 

Nh u = {\gl + ^gl) - 3Tr(Y ut n - Tr{Y^Y v ) 

N Hd = (lg 2 2 + ^g!)-3Tr(Y^Y d )-Tr(Y^Y e ) 



(15) 



iV s = -Tr{Yl R Y RR ) (16) 



where t — In// is the M5 1 scale) and 7 is the unit matrix. The RGEs for the gauge 
couplings are 

rln, 1 

bigf (17) 



dt 16tt 2 
where h = (f ,1,-3). 

The RGEs above are used to run the Yukawa matrices down from high energies 
(say the unification or string scale) down to the heaviest right-handed neutrino mass, 
M R z, which we assume is equal to the VEV of the S field. At this mass scale we 
perform a rotation of the right-handed neutrino fields to the basis in which M RR is 
diagonal, according to Eq.|5|, so that Y v is replaced by Y >v = Y U Q RR . In this basis 
one then runs the remaining RGEs (apart from Y RR ) down through the right-handed 
neutrino thresholds diag(M RZ , . . . , M R1 ), decoupling each right-handed neutrino at 



its mass threshold from the Y v contributions which appear on the right-hand side of 
the RGEs. To be explicit, we replace on the right-hand sides of the RGEs, 

K u - K U0 p ( 18 ) 

where 8 P = 6 (hi \i — hi Mr p ). For the (diagonal) Yrr on the right-hand side we replace 
it by 

Yrrpp yRRppOp- (19) 

In the next section we shall see that for the cases of interest it is not necessary to 
diagonalise Mrr in order to implement decoupling of the right-handed neutrinos, and 
one may remain in the basis defined by the U(l) family charges, since decoupling is 
facilitated by the simple structures of Mrr. 

Below the mass of the lightest right-handed neutrino one may use the RGE for k, 
the coefficient of the dimension 5 neutrino mass operator, 

~T + = -TTTaK 6 ^ + " ^Tr(Y u Y^) - (Y e Y^)n - k(Y^) t ] (20) 

at lo7r z 5 

In solving Eqj20|, it is convenient to diagonalise Y u and Y e at the scale of the lightest 
right-handed neutrino, using Eqs^, then make the approximation of keeping only the 
largest third family Yukawa eigenvalues. In the diagonal charged lepton basis, k must 
be transformed to k' given by 

«' = V gL kV} l (21) 

then, retaining the third family t and r Yukawa couplings only, the RGEs for the 
elements of k' are given by 

^ = [Q9l + pi - 6^ 2 - S i3 h 2 T - 5 3j h 2 T ] (22) 

Following from Eq.|2^ we see that the elements of m' LL (Mjn) = u|k'(Mri) at high 



energy are renormalised down to m' LL (m t ) = u|K'(m^) at low energy, ignoring the 
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running of v 2 , according to 



m' LL (m t ) = es /91 e 6/92 e 6It 



( m' LLU (M R1 ) m' LL12 (M R1 ) m' LL13 (M R1 )e^ 
m' LL21 {M m ) rn' LL22 (M R1 ) m , LL23 (M R1 )e~^ 
V m' LL31 (M m )e-^ m' LL32 (M R1 )e-^ m' LL33 (M R1 )e' 21 - 

(23) 



where 

I fin M R1 \ plnM m 



> = i5?jL. '«=iS?jL. 9lm (24) 



where f — t,r and i — 1, 2, 3. In running down to m t the charged lepton matrix 
will remain diagonal to good approximation, so that the low energy MNS matrix is 
simply given by 

Vmns — VJl- (25) 
where V' vL is the matrix which diagonalises m' LL (mt), 

m LL 9 ( m t) = Vl L m' LL {m t )V' v l = diag(m^ , m U2 ,m U3 ). (26) 

The mixing matrix Vmns and hence neutrino mixing angles, are the running 
quantities which can be computed at different energy scales. For example, the running 
of the neutrino mixing angle relevant to the atmospheric neutrino deficit, #23, can be 
understood from the evolution equation 

,33' 1 ^,22' 



rf „;Jofl _ o„:Jnfl n „:„2 a u2\ m LL + m LL 



16tt 2 - sin 2 26 23 = -2 sin 2 2^3 (1 - sin J 23 )(^ - K)-£ gr (27) 

at F tutt — 



This equation (see Babu et al in ref.Q) describes the evolution of the physical 23 
mixing angle assuming that we are in the diagonal charged lepton mass basis. 

3 Three right-handed neutrinos and SRHND 

We now specialise to three right-handed neutrinos and review the conditions for 
achieving SRHND and the LMA MSW solution g. The statement of SRHND is 
that, of the three right-handed neutrinos, one of them, N R3 , makes the dominant 
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contribution to the 23 block of tjill- This ensures that the 23 sub-determinant ap- 
proximately vanishes, and a 23 mass hierarchy therefore naturally results. 



We first write the neutrino Yukawa matrix in general as 

(a! a d \ 

Y v = b' b e (28) 
V d c f J 

There are now three distinct textures for the heavy Majorana neutrino matrix which 
maintain the isolation of the dominant right-handed neutrino Nrs, namely the diag- 
onal, democratic and off-diagonal textures introduced previously || [|. We consider 
each of them in turn. Note that, assuming SRHND, the contribution to the lepton 
23 and 13 mixing angles from the neutrino sector are approximately 

e d 



tan 9 



23 



tan# 



so that Super-Kamiokande and CHOOZ [THl imply 



(29) 



d < e « / 



(30) 



The condition on the 12 mixing angle such that it is relevant for the LMA MSW 
solution is discussed separately for each case below. 



3.1 Diagonal Texture 



M. 



RR 



( X' \ 

0X0 (31) 
V OF/ 

We may invert the heavy Majorana matrix and construct the light Majorana matrix 



using the see-saw mechanism, 



m LL 



I _i_ 9l- _l «f. 
' Y ~^ X X' 



V 



de _i_ ab i alb 1 

y„ K X> 
£ _l t. 4. k!l 

x ^ X' 



Y 



df I etc I a'c' \ 
y ~r X ~r X' I 
e£ , be , 6V 

y x ^ x' 

£ + + ^ 

y ~ x ~ x' 



The SRHND condition is 



e ef f xy x'y' 

rsj rs_/ ^>> 

y y y x' x' 



(32) 



(33) 



where x,y G a,b,c and x',y' G a',b',d. Q The 12 mixing angle determines whether 
we have the LMA MSW or SMA MSW solution, and this depends on the relative 
magnitude of the sub- dominant entries of m^, as discussed in |J. The condition for 
LMA MSW is M 



max 



ab ac a'b' a'd 
A' A 7 ' A 7 ' A 7 



'b' 2 be c 2 b' 2 b'c' c' 2 ' 
A' A' A' A 7 ' A 7 ' A 7 



max 



(34) 



3.2 Off-Diagonal Texture 



This is defined by: 

/ A \ 
M RR = AO (35) 
V OF/ 

The off-diagonal case is qualitatively different from the other two cases and gives 



/ d? _|_ 2aa' 



de I a'b i ab' df , a'c i ac' \ 

y~ l ~x~ t ~x y x x ' 



Y 



V 



y ' x 



A' 



SRHND is now defined by the conditions 

^ 2 e/ / 2 



e ej j xx 



where x G a, 6, c and x' G a', 6', c', The LMA MSW solution condition is 

(a'b ab' a'c ac'\ (W b'c be' cc'\ 

max — , — , — , — ~ max — , — , — , — 
\ X X X X \ X X X X 



(36) 



(37) 



(38) 



3.3 Democratic Texture 



The democratic case (assuming the Majorana masses in the upper block are of the 



same order but are not exactly equal) is defined by: 

/ A A 
M RR = A A 
V Y 



(39) 



2 The beauty of SRHND is that it automatically implies to„ 2 <C m V3 due to the approximately 
vanishing 23 subdeterminant, without the need for appeal to cancellations. To understand this 
simply drop the 1 / X terms and observe that the 23 subdeterminant vanishes which implies a massless 
eigenvalue which is a rather extreme case of a hierarchy! 
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The order of magnitude of uill is: 



/# + 0(f)+0(f) f + 0(f)+0(^) f + 0(f) + 0(^)\ 

£ + 0(§) + 0($) M + (|) + 0(f) 



(40) 



In this case the SRHND conditions are: 



e 2 e/ f 2 xy x'y' 

— ~ — ~ — ^> — ~ 41) 

Y Y Y X X K J 

where x,y G a,b,c and x', y' G a', b 1 , c'. The LMA MSW condition for a large 12 angle 
is 

fab ac a'V a'c'\ ( b 2 be c 2 b' 2 b'd c ,2 \ 

max [x ' x ' X' -x r max \ x' x ' x' x ' X' x (42) 



3.4 U(l) Family Symmetry 



Introducing a U(l) family symmetry | I4]| , ||15|| , [ 1^] , |p!7| provides a convenient way 



to organise the hierarchies within the various Yukawa matrices. For definiteness we 
shall focus on a particular class of model based on a single pseudo-anomalous U(l) 
gauged family symmetry |16| . We assume that the U(l) is broken by the equal VEVs 



of two singlets 8,8 which have vector-like charges ±1 fl6| . The U(l) breaking scale 



is set by < 8 >=< 8 > where the VEVs arise from a Green-Schwartz mechanism |18 
with computable Fayet-Illiopoulos D-term which determines these VEVs to be one or 
two orders of magnitude below Mjj. Additional exotic vector matter with mass My 
allows the Wolfenstein parameter [0 to be generated by the ratio [[16 

< 8 > <8> 



My My 



X w 0.22 (43) 



The idea is that at tree-level the U(l) family symmetry only permits third family 
Yukawa couplings (e.g. the top quark Yukawa coupling). Smaller Yukawa couplings 
are generated effectively from higher dimension non-renormalisable operators corre- 
sponding to insertions of 8 and 8 fields and hence to powers of the expansion parameter 

11 



in Eq.|43|, which we have identified with the Wolfenstein parameter. The number of 
powers of the expansion parameter is controlled by the U(l) charge of the particular 
operator. The fields relevant to neutrino masses Li, N£, H u , S are assigned U(l) 
charges l iy n p , h u = 0, a. From Eqs.|43|, the neutrino Yukawa couplings and Majorana 
mass terms may then be expanded in powers of the Wolfenstein parameter, 



M, 



RR 



/ ^|2tii+(t| \\ni+ri2+o-\ \\ni+ri3+a\ 

A |2n 2+ff | A |„ 2 +„ 3+CT | | <E> (44) 



The conditions which ensure that the third dominant neutrino is isolated require that 
the elements A l n i+ n 3+°1 ) A l n 2+™3+°-| sufficiently small. The diagonal, off-diagonal and 
democratic textures then emerge as approximate cases []8], p. The neutrino Yukawa 
matrix is explicitly 

^\h+ni\ \\h+n 2 \ A |Zi+n 3 | 

Y v ~ | Al /2+ni l Al' 2+n2 l Al /2+n3 l | (45) 

A lte+ni| ^\h+ri2\ ^\te+ri3\ 



which may be compared to the notation in Eq.28. The requirement of large 23 mixing 



and small 13 mixing expressed in Eq.^D| becomes 

\n 3 + h\ = \n 3 + l 3 \, \n 3 + k\ - \n 3 + l 3 \ = 1 or 2 (46) 

The remaining conditions for the U(l) charges depend on the specific heavy Majorana 
texture under consideration [f|. 

The charged lepton Yukawa matrix is given by 

(_\l'i+ei| \\h+e 2 \ \\h+e 3 \ \ 
A |i 2 +ei| A |i 2 +e 2 | A |/ 2 +e 3 | (47) 
\\h+£2\ \\h+e 3 \ J 

where e« are the U{\) charges of the charged lepton singlet fields. 

For the quarks we shall assume a common form for the textures of Y u and Y d 



I A 8 A 5 A 3 \ / A 4 A 3 A 3 

A 7 A 4 A 2 , Y d ~ A 3 A 2 A 2 | A n (48) 
V A 5 A 2 1 J V A 1 1 
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4 Renormalisation Group Analysis of SRHND 



In [0 we tabulated the simplest charges which satisfy all the conditions given above, 
and so provide a natural account of the atmospheric and solar neutrinos via the 
LMA MSW effect. In Table 1 we consider one example of each of the three cases, 
namely case A (Diagonal M RR ), case B (Off-diagonal M RR ) and case C (Democratic 
M RR ). The U(l) charges along with corresponding M RR , Y u , obtained from Eqs.pBI, 
[45], and other relevant parameters are outlined for each case, including the charged 
lepton charges, are also shown. Note that the zeroes in M RR appear after small angle 
rotations on the right-handed neutrino fields, which will not affect the perturbative 
expansion in powers of A in Y u . The order unity coefficients by, which are always 
present in U(l) models, are defined in this basis. 

In case A, before rotating to the diagonal charged lepton mass basis, we may 
easily estimate the order of the entries in in Eq.|32] from the matrices Y v and 
M RR in Table 1, and hence verify the SRHND conditions Eq.||, and the LMA MSW 
condition Eq.^4]. Using Eq.|32] we find 

/ A + A + A 1 + A + A 2 1 + A + A 2 \ 2 
m LL ~\ . i + A + A 3 i + A + A 3 (49) 

v ■ • i+A+A 3 ; <s> 

where the first entry in each element corresponds to the 1/Y contributions from N^ R , 
which clearly dominate the 23 block by a factor of A 2 , and the 12 and 13 elements 
by a factor of A. The 13 element is smaller than the elements in the 23 block by a 
factor of A leading to a 13 CHOOZ angle of this order. The subdominant entries in 
the 12,13,22,23 elements are the same order, leading to a large 12 angle suitable for 
LMA MSW. 

Similarly in case B, we can show that the SRHND condition Eq.[37] and LMA 
MSW condition Eq.RS are satisfied. We estimate the order of the entries in m^i in 
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Para- Case A 

meter (Diagonal Mrr 



Case B 

(Off-diagonal Mrr) 



Case C 

(Democratic Mrr) 



U(l) Z li2)3 = -1,1,1 
charges 711,2,3 = 1/2, 0, -1/2 
ei,2,3 = —3, —3, — 1, 
a= -1 



h,2,3 — 
™1,2,3 = 
e l,2,3 = 
(7 = 1 



-2,0,0 
-2,1,-1 
-2,-2,0, 



Z 1)2i3 = -1,1,1, 
"1,2,3 = 0, 0,-1/2 
ei,2,3 = —3, —3, —1, 
u — — 1 



/ anA 4 a i2 A 4 ai 3 A 2 

a 2 iA 2 a 22 A 2 a 2 3 
y a3iA 2 Q3 2 A 2 033 



( anA 4 ai 2 A 4 a i3 A 2 

a 2 iA 2 a 22 A 2 a 2 3 
\ a 3 iA 2 a 32 A 2 



/ anA 4 ai 2 A 4 CI13A 2 

a 2 iA 2 a 22 A 2 a 2 3 
V Q3iA 2 a 32 A 2 



ye 



«33 / 



a 3 3 / 



"V 



0.9 1.25 0.85 
1.2 4.2 1.25 
0.85 1.6 1.0 



0.9 1.25 0.85 
1.2 4.2 1.25 
0.85 1.6 1.0 




1.25 0.85 
4.2 1.25 
1.6 1.0 



&11A2 &i 2 A 613X2" 

3 1 

6 2 iAa 6 22 A 6 2 3A2 

3 1 

y &31A2 fo 32 A &33A2 



yv 



( 6 n A 4 6 12 A 6 13 A 3 
6 2 iA 2 6 22 A 6 2 3A 
\ 631 A 2 6 32 A 633 A 



611A b 19 \ &13A 



n2 



13' 



o 2 iA b??, A 6^3 A 



'22 



23' 



&31A 6 32 A 633 A 



33' 



0.5 0.85 1.0 

1.0 1.3 0.4 

1.1 0.4 1.5 




/ 0.5 0.85 1.0 
1.0 1.3 0.5 
v 1.1 0.5 1.5 



M RR 



/ 1 







A 
A 2 




/ A c i2 A 
c i2 A A 
V A 2 



< E > 



22.75 x VQ^GeV 
1.10 x 10 14 Ge\/ 



2.20 x lO^el/ 
4.84 x W A GeV 



38.20 x l& A GeV 
1.85 x 10 14 Ge\/ 



M 



Ri 



'91,2 



0.048918,0.0751076 
0.133729,0.118012 
4.1884.10" 4 , 1.794.10- 8 



0.046789,0.072753 
0.130667,0.111514 
4.0474.10" 4 , 1.711. 10- 8 



0.0498926,0.07625 
0.135464,0.120819 
4.278.10" 4 , 1.829.10- 



Table 1 : Textures of the Yukawa couplings of Dirac neutrino mass, right-handed Majorana neutrino 
mass, and also other relevant parameters (as defined in the text) required for the numerical estimation 
of left-handed Majorana neutrino masses at low energies through see-saw mechanism in tables 2,3 
and 4. < S > is taken as a free parameter and Mr x is the lowest threshold scale in Mrr- We take 
C12 « 0.9. 
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Eq.|36| from the matrices Y v and Mrr in Table 1 as, 



A 5 + 2A 5 A 3 + A 5 + A 3 A 3 + A 5 + A 3 
A + 2A 3 A + A 3 + A 3 
A + 2A 3 



) 




(50) 



where the first entry in each element corresponds to the 1/Y contributions coming 
from the right-handed neutrino N^r, and clearly dominates the 23 block by a factor 
of A 2 . In this case it does not dominate the other elements outside the 23 block. The 
13 element from N 3 r is suppressed relative to the 23 block elements by a factor of A 2 , 
leading to a CHOOZ angle of this order. The subdominant entries in the 12,13,22,23 
elements are again the same order, leading to a large 12 angle suitable for LMA MSW. 

Finally in case C, we verify that the SRHND condition Eq.|| and LMA MSW 
condition Eq.|42| are satisfied, by constructing the entries in in Eq.[40] from the 
matrices Y u and Mrr in Table 1, 



where the first entry in each element corresponds to the 1/Y contributions coming 
from the dominant right-handed neutrino, which dominates in the 23 block by a 
factor of A 2 and in the 12 and 13 elements by a factor of A. The 13 element is smaller 
than the 23 elements by a factor of A leading to a CHOOZ angle of this order. The 
subdominant 12,13,22,23 elements are all the same order leading to LMA MSW. 

We now turn to a numerical treatment of these three cases. In integrating the 
RGEs down to low energies the heavy thresholds for the diagonal texture in Eq.|3~T| 
are dealt with exactly as described in the previous section, except that no rotation 
is required to get to the diagonal Mrr basis since we already begin with that form. 
For the democratic texture in Eq.|39] there are essentially only two mass thresholds to 
consider X and Y since Nri and Nr2 are approximately degenerate with mass X and 
we can ignore any small mass difference between them to leading order. Similarly for 



/ A+A+A 1+A+A 1+A+A 
i+A+A i+A+A 
V • • i+A+A 




(51) 
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the off-diagonal texture in Eq.^ we also have only two mass thresholds to consider 
X and Y since Nm and Nr2 are now exactly degenerate with mass X. Thus for both 
democratic and off-diagonal textures we can replace on the right-hand sides of the 
RGEs by 

12-^0!, Y»^Y»0 2 , Y£^Y»9 3 (52) 

where 6*1,2 = #(ln/i — lnX), 6*3 = 6*(ln/i — In V). We replace Yrr on the right-hand 
side of the RGEs by 

YRRij —> Y RRi j9i9j. (53) 

Tables 2-4 give the numerical values of quantities at three different energy scales: 
the GUT scale Mjj = 2.0 x 10 16 GeV, the lightest right-handed neutrino mass 
~ 10 14 GeV, and the low energy scale m t = 175 GeV. To begin with the left- 
handed Majorana neutrino mass matrix 0^ll(^x) arising from the sum of dominant 
and subdominant contributions from Nr3 and Nri,Nr2 respectively, are numerically 
computed at the energy scale /1 = Mu = 2.0 x 10 16 GeV, as seen in tables 2,3,4 for 
the three cases (A,B,C) for large tan/3. The corresponding MNS mixing matrices 
which in turn give S so i and S at along with neutrino masses [rn v \,m V 2,rn v z] are also 
estimated. Next, as already discussed in section 2, we calculate the radiative correc- 
tions to m' LL and Vmns through the running of the Yukawa couplings Y u , Y d , Y u , 
Y e , Yrr, and the gauge couplings g±, g% from high energy scale Mjj through the 
heavy neutrino threshold region at successive steps down to the lightest right-handed 
neutrino mass (MrJ. The corresponding m' LL (M Rl ) (in the diagonal charged lepton 
basis) and Vmtvs^MrJ are again estimated as shown in tables 2,3 and 4 for three 
cases. Further, the radiative corrections to m! LL and Vmns from the lightest right- 
handed neutrino down to low energies ( say top-quark mass scale m t ) are taken in 

the usual way through the running of the coefficient of the dimension 5 operator n' 

3 Strictly speaking the see-saw mechanism does not operate at energy scales higher than the 
right-handed neutrino mass scales, but such results do provide a meaningful measure of the effects 
of radiative corrections from the GUT scale down to low energies. 
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Scale fi = M v = 2.0 x 10 16 GeV 



Scale /I = M v = 2.0 x 10 Lb GeV 



Y e = 

2.114.10- 3 
5.808.10" 2 
4.114.10- 2 



2.928.10- 3 

0.20328 

7.744.10- 2 



4.114.10- 2 \ 
1.250 
1.000 



Yliaa = ^(4.36.10- 4 , 6.61. KT 2 , 1.62) 



diag 
™>LL = 

( 4.192.10- 3 
7.724. 10- 3 

-2 



7.724.10- 3 1.534.10- 2 
5.345.10- 2 6.894.10- 2 
V 1.534.10- 2 6.894.10- 2 9.965.10~ 2 
m f^ = diag(2A6 x 10" 4 , 0.00582, 0.151) 



V eL = 

f 0.999 0.002 -0.044 \ 

0.036 -0.621 0.783 
V 0.026 0.784 0.620 / 
m e /m^,m^/m T = 0.0066,0.041 

VmNS = VeLV uL = 

0.797 -0.593 0.115 
0.432 0.692 0.579 
-0.422 -0.412 0.808 
S sol = 0.9173, = 0.8965 



Scale /i = M Rl = 1.10 x 10 l4 GeV 



Scale n = M Rl = 1.10 x 10 1A GeV 



Y e = 

1.411.10 3 1.358.10 
4.234.10- 2 0.15298 
2.890.10- 2 4.766.10- 
F^ = ^(3.74.10- 4 ,5 

m LL = 

3.827.10- 3 6.800.10- 
6.800.10- 3 4.848.10- 
1.278.10" 2 5.755.10" 



m '*w = dio ff (2.294.10- 



3 2.621. 10- 2 
0.90295 

2 0.71419 
72. 10- 2 , 1.16) 

3 1.278.10 -2 \ 
2 5.755.10- 2 

2 7.715.10- 2 / 
,0.0054,0.1238) 



V eL = 

( 0.996 0.004 -0.042 
0.035 -0.615 0.787 

V 0.023 0.788 0.615 
m e /m^,m^/m T = 0.0065,0.049 
Vmns = V eL V} L = 

I 0.790 -0.602 0.118 
0.421 0.671 0.606 

V -0.446 -0.432 0.784 
S sol = 0.92991, S at = 0.9366 



Scale fj, — rrtt — 175GeV 



Scale n = m t = 175GeV 



Yliaa = ^(2.89.10" 4 ,4.41.10- 2 , 0.627) 



m 



diag 

IV 

LL 



( 2.855.10" 3 5.073.10 
5.073.10- 3 3.616.10 
V 8.474.10- 3 3.816.10 
m'fi* = diag(\m x 10 



3 8.474. 10- 3 
2 3.816.10- 2 
" 2 4.545.10- 2 
4 , 0.003848, 0.08046) 



0.0065,0.0703 



Vmns = VyL = 
( 0.768 -0.628 0.124 
0.412 0.633 0.656 
V -0.490 -0.453 0.745 
S sol = 0.9606, = 0.9840 



Table 2: Results for case A. Left-handed Majorana neutrino mass matrix in the diagonal charged 
lcpton basis and the MNS mixing matrix at different energy scales M[/,M^ 1 ,m t for the case of 
diagonal Mrr. Neutrino masses are expressed in eV. 
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Scale n = M v = 2.0 x 10 w GeV 



Scale = Mu = 2.0 x 10 m GeV 



Y e = 

2.114.10" 
5.808.10- 
4.114.10" 



2.928.10" 

0.20328 

7.744.10" 



4.114.10- 2 \ 
1.250 
1.000 



Ydiaa = dtag(4.36.W-\ 6.61.10" 2 , 1.62) 



-0.044 \ 
0.783 
0.620 

0.0066,0.041 



/ 4.457.10- 5 

-9.965.10~ 4 
V 5.172.10- 3 
m LL = diag{- 



-9.965.10" 4 5.172.10- 3 
3.116.10- 2 3.784.10~ 2 
3.784.10~ 2 5.983.10- 2 
-0.21 x 10- 2 , 0.00701, 0.086) 



V eL = 

f 0.999 0.002 

0.036 -0.621 
V 0.026 0.784 

m e /m^m^/m T = 

VmNS = VeLV uL = 

0.886 -0.462 0.042 
0.363 0.740 0.566 
-0.296 -0.485 0.823 
S sol = 0.6731, = 0.8918 



/ 



Scale (j, — M Rl — 4.84 x 10 13 GW 



Scale n = M Rl = 4.84 x 10 u GeV 



Y e 



ye 
diag 



1.358.10- 3 1.166.10 3 2.527.10~ 2 
4.072.10- 2 0.148 0.868 
2.803.10- 2 4.548.10- 2 0.6940 

diag(3.Q7.10~\ 5.64.10~ 2 , 1.12) 



-0.043 \ 
0.784 
0.619 J 
0.0065,0.0503 



m 



m 



LL — 

-1.733.10~ 5 
— 1.113.10 3 
4.481. 10- 3 
ll = diag{- 



-1. 113.10- 3 

2.972.10- 2 

3.333.10- 2 



4.481.10" 
3.333.10" 
4.907.10" 



■2.035.10~ 3 , 0.00659, 0.0/ 



V eL = 

/ 0.999 0.005 

0.037 -0.620 

V 0.023 0.785 

VmNS — VelYvL = 

0.882 -0.469 0.039 \ 
0.355 0.718 0.598 
-0.309 -0.514 0.800 / 
S sol = 0.6875, = 0.9201 



Scale ji — rrit— 175GeV 



Scale fx — mt— 175GeV 



= diag(2.85.1Q-\ 4.39. 10~ 2 , 0.619) 



diag 

™'L = 

( -1.295.10~ 5 
-8.317.10~ 4 
V 2.995.10- 3 
mf L a9 = dzag{- 



-8.317.10~ 4 2.995.10 
2.221. 10~ 2 2.228.10 
2.228.10" 2 2.934.10 
1.486 x 10~ 3 , 0.00462,0 



-2 
-2 



0484) 



m e /m^mJm T = 0.0065,0.0709 

VmNS = VuL = 

0.876 0.480 0.036 
0.346 -0.680 0.647 
-0.335 0.555 0.762 
S sol = 0.7102, S at = 0.9738 



Table 3: Results for case B. Left-handed Majorana neutrino mass matrix in the diagonal charged 
lcpton basis and the MNS mixing matrix at different energy scales Mjj, Mn x ,mt for the case of 
off-diagonal Mrr. Neutrino masses are expressed in eV. 
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Scale n = Mu = 2.0 x 10 m GeV 



Scale /I = Mu = 2.0 x 10 Lb GeV 



Y e = 

2.114.10- 3 
5.808.10- 2 
4.114.10- 2 



2.928.10- 3 

0.20328 

7.744.10- 2 



4.114.10- 2 \ 
1.250 

1.000 / 



YL g = ^(4.36.10- 4 ,6.61.10- 2 ,1.62) 
m'E L = 

f 2.751.10- 3 2.878.10- 3 8.715.10- 3 
2.878.10- 3 3.070.10- 2 4.222.10- 2 
V 8.715.10- 3 4.222.10-2 6.736.10- 2 
m fw = diag(2.9Q x 1Q- 5 , 0.00488, 0.096) 



V eL = 

( 0.999 0.002 -0.044 \ 

0.036 -0.621 0.783 
V 0.026 0.784 0.620 / 
m e /m^,m^/m T = 0.0066,0.041 

VmNS = VeLV uL = 

0.780 -0.619 0.095 
0.480 0.688 0.544 
-0.402 -0.379 0.834 
S sol = 0.9488, = 0.8378 



Scale n = M Rl = 1.85 x 10 l4 GeV 



Scale fi = M Rl = 1.85 x 10 l4 GeV 



ye 



1.391.10 3 1.274.10 
4.293.10-2 0.1550 
2.935.10-2 4.864.10 



ye 

diag 



3 2.572.10- 2 

0.9159 
2 0.7248 
diag(3. 77. lQ-\ 5.77.10- 2 , 1.18) 



m 



LL 



2.433.10" 3 2.426.10" 

2.426.10- 3 2.753.10- 

7.009.10- 3 3.471.10- 

/diag 



3 7.009.10- 3 \ 
2 3.471.10-2 
2 5.094.10- 2 j 



m'l™9 = diag(2.Q7A.10~ 5 , 0.00433, 0.07647) 



V eL = 

( 0.999 0.005 -0.041 
0.035 -0.616 0.787 

V 0.022 0.788 0.615 
m e /m^,m^/m T = 0.0065,0.049 
Vmns = V eL Vl L = 

f 0.771 -0.630 0.096 \ 
0.472 0.665 0.578 

V -0.428 -0.401 0.810 / 
S sol = 0.9604, S,a = 0.8945 



Scale fj, — rrit — 175GeV 



Scale (j, — m t — 175GeV 



Yf 



diag -= ^(2.89.10" 4 ,4.41.10- 2 , 0.626) 
m'E L = 

( 1.811. 10- 3 1.805.10 
1.805.10- 3 2.049.10 
V 4.623.10" 3 2.289.10- 2 
m'%* = diag(1.90 x 10 



3 4.623.10- 3 
2 2.289.10-2 
2.978.10-2 
5 , 0.00309, 0.04897) 



m e /m^,m^/m T = 0.0065,0.0703 

Vmns — V vL = 
/ 0.752 -0.652 0.099 
0.461 0.628 0.627 
V -0.471 -0.426 0.772 
S sol = 0.9798, Sat = 0.9579 



Table 4: Results for case C. Left-handed Majorana neutrino mass matrix in the diagonal charged 
lcpton basis and the MNS mixing matrix at different energy scales M[/,M^ 1 ,m t for the case of 
democratic Mrr. Neutrino masses are expressed in eV. 
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in the diagonal charged lepton basis. These low energy results in tables 2,3,4 can be 



compared with observational data, and S so i, S at defined in Eqs. |T3| , |T4] and the neutrino 
masses given as the diagonal elements of m^™ 9 are seen to be in good agreement with 
atmospheric data and the LMA MSW solutions. The CHOOZ constraint is also satis- 
fied in all cases, with the 13 element of Vmns being larger for cases A and C than for 
case B, as expected from the analaytical estimates above. The charged lepton mass 
ratios are rrie/m^ pa 0.0066 and m^/rrir pa 0.07 at low energy scale, compared to the 
experimental values 0.005 and 0.06, respectively. 

It is instructive to examine the RG evolution of the neutrino masses and mixing 
angles. The numerical results in Tables 1-4 show that neutrino masses m U2 and m U3 
are decreasing from high energy scale Mjj to low energy scale m t by about 40%, but 
the ratios m U2 /m V3 increase by about (24%, 17%, 55%) corresponding to low energy 
ratios of (0.05, 0.10, 0.08) for cases (A,B,C), respectively. The atmospheric mixing 
quantity S a t increases by about 10% in the three cases, and approaches maximal 
mixing in all cases. This can be understood from Eq.^7| which shows that in the 
diagonal charged lepton basis, the sign of the RGE evolution is determined by the 
sign of m^ L — m 2 ^' L) which in the examples in Tables 2,3,4 is always positive. This 
means that the overall sign of the RGE is negative, which implies an increasing S at 
as the energy scale is reduced. Note that the largest contribution to the atmospheric 
mixing angle arises from the charged lepton sector, for this choice of parameters, 
as is clear from examining V e i in Tables 2,3,4. The solar mixing quantity S so i also 
increases, but the increases of 1 — 5% are mild in comparison to S a t- In cases A,C 
Ssoi approaches maximal mixing, while in case B its low energy value is about 0.71, 
which are well near the best fit sir?1Q\2 ~ 0.76. 

In Figures 1-3 we display the RG running of the mixing angles S so i and S at and 
the ratio of neutrino masses m v ijm v z as a function of t = ln/i, for the three cases 
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A,B,C. We prefer to show the variation in the neutrino mass ratio, rather than the 
absolute values of the two neutrino mass eigenvalues since they will depend to some 
extent on the running vacuum expectation value, although the effect on the present 
analysis is not large since we assume high tan (3. We also do not consider the ratio 
m v i/m V 2 since this is not experimentally measurable in the hierarchical case. 

In Figure 1(a) we show the variations of S so i and S at with energy scale for case 

A. As already stated, the atmospheric mixing angle runs more rapidly than the solar 
mixing angle, and although S at starts out smaller then S so i at Mjj, it quickly grows 
larger. Note the effect of the heavy right-handed neutrino mass thresholds which 
change the slope of the curves, which grow steeply over the heavy threshold region. 
In Figure 1(b) we give the variation of the neutrino mass ratio m v ifm v % with t — In// 
for case A. Here the effects of the three heavy right-handed neutrino thresholds is 
clearly seen, with again a steep rise in this mass ratio over the heavy threshold region. 

In Figure 2(a) we show the variations of S so i and S at with energy scale for case 

B. Here the atmospheric angle starts out larger than the solar angle, but still grows 
more rapidly. In Figure 2(b) we give the variation of the neutrino mass ratio m u2 /m u3 
with t — In// for case B. The qualitative shape of this curve is similar to Figure 1(b), 
but there are only two heavy neutrino mass thresholds in this case, and also the mass 
ratio is larger throughout. 

In Figure 3(a) we show the variations of S so i and S at with energy scale for case B. 
Here the solar angle starts out much larger than the atmospheric angle, and as before 
the atmospheric angle grows more rapidly and approaches the solar angle. In Figure 
3(b) we give the variation of the neutrino mass ratio m v2 /m„3 with t — In// for case 
B. Because of the choice of c± 2 in Table 1, the two heavy right-handed neutrino mass 
thresholds are very close together in this case. 
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Fig. 1(a) Case A (Diagonal Mrr): Variation of S so i and S a t with energy scale t 
In /i, which are represented by solid-line and dotted-line respectively. 
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Fig. 1(b) Case A (Diagonal Mr#): Variation of the neutrino mass ratio m„ 2 /m„3 
with t = In /i. 
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Fig. 2(a) Case B (Off-diagonal Mrr): Variation of S so i and S a t with t = ln/x, which 
are represented by solid-line and dotted-line respectively. 
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Fig. 2(b) Case B (Off-diagonal Mrr)- Variation of the neutrino mass ratio m„ 2 /m„3 
with t = In /i. 
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Fig. 3(a) Case C (Democratic Mrr): Variation of S so i and S a t with energy scale 
t = ln/x, which are represented by solid-line and dotted-line respectively. 
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Fig. 3(b) Case C (Democratic Mrr): Variation of the neutrino mass ratio m^/m^ 
with t — ln/x. 
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It is important to emphasise that more than 50% of the radiative corrections to 
all the physical quantities m U2 /m U3 , S at , S so i arises from the RG evolution over the 
range from the GUT scale Mjj = 2.0 x 10 16 GeV through the heavy right-handed 
neutrino threshold region down to the lightest right-handed neutrino mass scale 
Mrj ~ 10 14 GeV. Therefore even though this range covers only two orders of magni- 
tude in energy, it is just as important as the RG evolution effects from M Rl ~ 10 14 GeV 
down to low energies which covers over 12 orders of magnitude in energy, and is the 
region commonly considered in the literature. 

5 Conclusion 

We have studied the effects of radiative corrections on neutrino masses and mixing 
angles, in evolving the theory defined at the GUT scale down to low energies. Our 
approach is to run down all the Yukawa matrices from the GUT scale down through 
the heavy right-handed neutrino mass thesholds to low energy, replacing the neutrino 
Yukawa matrices by the dimension 5 neutrino mass operator at the lowest right- 
handed neutrino mass threshold. We have found that in the realistic cases considered, 
the atmospheric and solar neutrino mixing parameters receive radiative corrections 
of order 10% and 5%, respectively, while the ratios of the neutrino masses change by 
up about 50% in going from My to low energy. Importantly, more than 50% of the 
overall radiative corrections arises from running through the threshold region, even 
though it only accounts for two orders of magnitude in energy. Thus many of the 
existing analyses in the literature which ignore this threshold region could endanger 
a significant error. 

We have considered a realistic class of models known as SRHND [§, |^| , in which the 
contribution to the 23 block of the light effective Majorana matrix is dominated 
by a single right-handed neutrino. The small neutrino mass hierarchy m U2 /m U3 ~ 0.1 
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originates from the smallness of the 23 sub determinant of tull (which is zero in 
the limit that a single right-handed neutrino is the only contribution). We have fo- 
cussed on cases with two large mixing angles #23 and #i 2 , with the CHOOZ angle #13 
being small - the so-called bi-maximal mixing scenario. The couplings of the domi- 
nant right-hand neutrino controls the 23 and 13 mixing angles, and the subdominant 
right-handed neutrino couplings control the 12 angle. A U(l) family symmetry is 
used to generate a controlled expansion of all the Yukawa couplings in powers of the 
Wolfenstein parameter A, and with a suitable choice of U(l) charges SRHND may be 
achieved with subdominant right-handed neutrino contributions being of the correct 
order of magnitude to generate the desired neutrino mixing angles and mass hierarchy. 

The U(l) charges also control the charged lepton Yukawa matrix, resulting in 
large contributions to the physical lepton mixing angles from the charged lepton 
sector. In general one would expect the two contributions to the 23 mixing angle 
(from the charged leptons and the neutrinos) not to cancel exactly due to order unity 
coefficients in the Yukawa matrices which are not predicted, and our numerical results 
include these effects. Similarly it would be surprising if the physical mixing angles 
turned out to be exactly maximal at the GUT scale. Therefore we have considered 
examples in which the 23 and 12 angles start out large but not maximal, and we 
have also assumed large tan/5, where the r Yukawa coupling is large and the effect 
of radiative corrections is maximised. Although the effect of radiative corrections on 
the mixing angles is always < 10%, showing that the models are quite stable, we 
have shown that the effects may play an important role in driving the initially large 
(but not maximal) mixing angle towards its maximal value. This is true both of the 
atmospheric mixing angle and the solar mixing angle, although in the latter case the 
effects are milder, which helps to explain why the atmospheric angle is larger than 
the solar angle. In principle, a different choice of parameters could have caused the 
neutrino angles to have grown smaller at low energies. However it is significant that 
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for the cases considered both mixing angles become magnified showing that the low 
energy approximate bi-maximal scenario could partly result from radiative corrections 
in SRHND models. 
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